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As  waves  propagate  into  shallower  waters,  their  profile  changes  fi'om  a  Gaussian 
distribution  to  an  increasingly  non-Gaussian  profile.  The  non-Gaussian  properties  of 
waves  in  finite  water  depths  are  analyzed  through  the  use  of  a  probability  density  fimction 
defined  in  closed  form  and  that  can  be  solved  by  the  application  of  wave  displacement 
data.  Initially,  the  applicability  of  the  probability  density  fimction  is  verified  using  wave 
data  records  obtained  from  the  Coastal  Engineering  Research  Center  (CERC)  at  Duck, 
North  Carolina.  A  broad  range  of  sea  conditions  and  water  depths  are  represented  to 
ensure  the  probability  density  fimction  correlates  well  with  actual  histogram  data  over 
varied  conditions.  With  the  applicability  verified,  the  probability  density  fimction  is  used 
to  define  criteria  for  the  boundary  where  the  wave  field  can  no  longer  be  considered 
Gaussian.  The  ability  to  determine  when  a  shoaling  wave  profile  can  no  longer  be 
considered  Gaussian  is  of  paramount  inqrortance  in  order  to  ensure  proper  wave  theories 
are  being  applied.  This  boundary  between  Gaussian  and  non-Gaussian  wave  profiles  will 
be  determined  in  terms  of  significant  wave  height  and  water  depth. 
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Presently,  wave  data  records  are  required  to  define  the  probability  density  fiinction 
used  in  this  study.  Trend  anafyses  were  conducted  on  the  parameters  of  the  probability 
density  fimction  to  determine  if  shiqilifications  can  be  made  to  its  form  and  to  see  if  it  is 
possible  to  define  the  probability  density  fimction  in  terms  of  the  local  conditions,  sea 
severity  and  water  depth,  without  the  need  for  wave  data  records. 


CHAPTER  1 
INTRODUCTION 

As  wind  generated  waves  transition  jfrom  deep  water  to  shallow  water,  a  significant 
change  is  evident  in  the  wave  profile.  As  an  exanqple.  Figure  1-1  shows  portions  of  wave 
records  obtained  by  the  Coastal  Engineering  Research  Center,  (CERC),  during  the 
ARSLOE  Project.  Figure  l-l(a)  is  the  wave  profile  recorded  at  a  distance  of  60  meters 
of^ore  in  a  shallow  intermediate  water  depth  condition.  Figure  1-1  (b)  is  the  wave 
profile  taken  at  the  same  time  but  at  a  distance  of  12  km  offihore  in  a  deep  water 
condition.  Deep  water  wave  profiles  in  the  sea  characteristically  reflect  a  Gaussian  (linear) 
random  process.  The  height  of  the  wave  crests  are  probabilistically  equal  to  the  height  of 
the  wave  troughs.  As  waves  move  from  deep  to  shallow  water,  the  wave  profile 
transforms  from  a  Gaussian  random  process  to  a  non-Gaussian  (nonlinear)  random 
process  in  which  the  wave  profile  shows  a  definite  excess  of  high  crests  and  shallow 
troughs. 


a)  Gage  615,  Severe  Sea  Condition  b)  Gage  625,  Mild  Sea  Condition 


Figure  1.1  Wave  Profile  Examples 
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TMs  remarkable  transformation  in  the  wave  profile,  as  an  irregular  sea  transitions 
fi'om  deep  to  shallow  water,  has  also  been  evidenced  in  laboratory  experiments  (Doering 
and  Donelan  1993).  The  transformation  of  wave  profiles  can  be  attributed  to  the  non¬ 
linear  wave-wave  interactions  augmented  primarily  by  decreasing  water  depth.  The 
“degree”  of  non-linearity  a  wave  field  has  at  a  particular  site  of  interest  in  the  coastal 
region  appears  to  depend  largely  on  sea  severity  and  water  depth. 

Although  there  has  been  awareness  that  Gaussian-based  solutions  do  not  always 
represent  accurate  resuhs  in  the  coastal  region,  clarification  of  wave  characteristics  fi'om 
deep  to  shallow  water  has  been  a  subject  of  considerable  interest  in  ocean  and  coastal 
engineering.  There  have  been  few  studies  which  have  offered  non-Gaussian  solutions 
describing  the  probability  density  fimction  of  waves  in  coastal  regions.  A  commonly 
known  study  was  by  Longuet-Higgins  (1963),  which  statistically  predicts  coastal  wave 
properties  using  a  power  series  ejq)ansion.  This  probabihty  distribution  is  based  essentially 
on  the  Gram-Charlier  series  distribution.  A  probability  density  fimction  developed  by  Ochi 
and  Ahn  (1994a)  appears  to  be  applicable  throughout  the  coastal  region  fiom  deep  to 
shallow  water  and  is  defined  in  closed  form  Equation  (1). 

fi;x)  =  -==L — exp  -  — — ^ — ^(1  -  -  exT^{-yax)f  -  yax 

\27t<7.  L  2{yacT,) 

.  L28^r:x>0 
where,  y  =  \ 

^  3.00/or:x<0 

This  is  the  probabihty  density  function  that  will  be  used  as  the  basis  to  analyze 
coastal  wave  profiles  for  this  study. 


)ose  of  the  Stud\ 


The  overall  objective  of  this  study  involves  two  components.  The  first  component 
will  be  to  vaUdate  the  apphcabihty  of  the  probabihty  density  function  developed  for  waves 
in  finite  water  depth  (Ochi  and  Ahn  1994a)  and  to  analyze  the  transition  point  vsfiere  wave 
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profiles  can  no  longer  be  considered  Gaussian.  In  order  to  acconq)lish  the  first  objective  of 

this  study,  a  statistical  anafysis  was  carried  out  on  wave  time  histories  obtained  during  a 

storm  by  the  Coastal  Engineering  Research  Center  (CERC)  during  the  Atlantic  Ocean 

Remote  Sensing  Land-Ocean  Experiment  (ARSLOE)  Project  in  1980.  To  interpret  the 

transformation  of  waves  fi’om  a  Gaussian  to  a  non-Gaussian  random  process,  wave 

records  at  seven  locations  were  taken  where  the  water  depths  ranged  fi'om  two  meters  to 

twenty-one  meters.  In  total,  approximately  1000  wave  records  were  analyzed  with  each 

wave  record  consisting  of 4800  water  displacement  measurements.  The  individual  wave 

» 

displacement  measurements  were  recorded  at  a  fi’equency  of  4  Hertz  in  twenty  minute 
intervals. 

Upon  confirmation  of  the  prohabihty  density  fimctions’  apphcability,  the  histogram 
data  will  be  con^ared  closely  with  the  corresponding  normal  distribution  to  determine  the 
conditions  imder  which  the  wave  profile  can  no  longer  be  considered  Gaussian.  Li 
particular,  general  relationships  to  sea  severity  and  water  depth  will  be  determined.  A 
criteria  will  be  established  which  will  define  the  limit  to  when  the  Gaussian  wave  profile 
becomes  non-Gaussian.  This  limit  will  be  based  on  the  nonlinearity  parameters  analyzed 
fi'om  the  wave  profile  data.  The  second  objective  of  the  study  will  then  be  to  analyze  the 
parameters  of  the  probabihty  density  fimction.  The  three  parameters,  a,  and  a*  will  be 
statistically  analyzed  in  both  dimensional  and  non-dimensional  form  Trends  based  on  the 
analysis  will  be  described  and  relationships  linking  these  parameters  of  the  probabihty 
density  fimction  to  site  conditions  wiU  be  developed. 

This  study  consists  of  six  chapters.  Chapter  2  provides  a  general  review  on  the 
development  of  probabihty  density  fimctions  apphcable  to  coastal  regions.  Chapter  3 
describes  the  ARSLOE  Project  of  1980  and  addresses  the  statistical  approach  taken  in  this 
study.  The  data  used  as  the  basis  for  this  study,  were  obtained  fi'om  the  ARSLOE  Project. 
This  chapter  whl  also  describe  the  local  field  conditions  at  CERC  diuing  the  time  the  data 
were  cohected.  Chapter  4  addresses  the  probabihty  density  fimction  apphcabihty  and 
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analysis.  This  chapter  has  two  distinct  sections.  The  first  section  reviews  the  finite  water 
depth  probability  density  fimction  (Ochi  and  Ahn  1994b)  and  verifies  the  applicability  of 
the  probability  density  fimction.  The  second  section  discusses  sea  severity  and  water 
depth  and  the  transition  between  Gaussian  and  non-Gaussian  wave  profiles.  Chapter  5 
relates  the  non-dimensionalized  parameters  of  the  probability  density  fimction,  aa,  |ii»/a 
and  0*1(5,  to  sea  severity,  (a),  and  water  depth  in  both  dimensional  and  non-dimensional 
form  The  parameters  will  also  be  cort^ared  to  each  other  and  relationships  determined. 
This  chapter  describes  trends  evident  based  on  statistical  analysis.  These  trends  are 
qualitative  and  are  intended  to  give  a  better  understanding  of  the  interaction  the  wave 
profile  has  with  its  environment  as  it  propagates  in  shallow  waters.  Finally,  Chapter  6 
provides  the  summary  and  conclusion  for  this  study. 


CHAPTER  2 
LITERATURE  REVIEW 


As  stated  in  the  Litroduction,  the  profile  of  waves  in  finite  water  depths  is  quite 
different  fi:oni  that  observed  in  deep  water.  From  the  stochastic  view  point,  waves  in  finite 
water  depth  are  categorized  as  a  non-Gaussian  random  process  for  \\Mch  a  precise 
description  of  the  process  does  not  exist.  This  is  in  contrast  to  the  Gaussian  random 
process  which  can  be  defined  by  a  rigorous  mathematical  formulation.  Although 
considerable  attention  has  been  given  to  Gaussian  waves  in  deep  water,  little  information 
is  available  on  the  non-Gaussian  characteristics  of  shallow  water  waves. 

One  of  the  most  commonly  known  presentations  describing  non-Gaussian  waves  is 
the  Gram-Charlier  series  probabihty  density  fimction.  This  density  fimction  is  given  in  a 
series  consisting  of  an  infinite  number  of  terms.  Longuet-Higgins  (1963)  theoretically 
showed  that  the  statistical  properties  of  nonlinear  waves  can  be  presented  by  the  Gram- 
CharMer  series  distribution.  Bitner  (1980)  found  reasonable  agreement  in  a  coniparison 
between  the  Gram-CharUer  series  distribution  and  the  histogram  constructed  fi’om  wave 
data  obtained  in  a  relatively  shallow  water  area. 

On  the  other  hand,  Ochi  and  Wang  (1984)  analyzed  more  than  500  wave  records 
obtained  by  the  Coastal  Engineering  Research  Center  (CERC)  at  Duck,  North  Carolina 
during  the  ARSLOE  Project  (1980).  They  fmmd  that  while  the  statistical  properties  of 
waves  at  any  water  depth  can  be  well  represented  by  the  Gram-Charher  series  distribution, 
there  is  a  drawback  to  applying  this  density  fimction  in  practice.  The  probability  density 
fimction  becomes  negative  for  large  negative  wave  displacements  and  this  negative  portion 
does  not  vanish  with  increasing  niunber  of  terms  in  the  series  distribution.  Because  of  this 
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difficulty,  the  transition  property  of  wave  characteristics  from  gaussian  to  non-gaussian 
distribution  can  not  be  well  clarified.  They  derived,  however,  a  limiting  sea  severity  for  a 
specified  water  depth  below  A^diich  wind-generated  coastal  waves  may  be  considered  to  be 
a  Gaussian  random  process. 

Another  approach  for  describing  the  surfice  elevation  of  non-Gaussian  random 

waves  is  to  assume  that  individual  waves  can  be  expressed  as  a  Stokes  wave.  There 

appears  to  be  some  reservation  to  applying  this  approach  for  representing  the  statistical 

distribution  of  the  wave  profile.  This  method  presupposes  a  preliminary  form  of  the  wave 

» 

profile  before  establishing  its  probability  density.  This  supposition  may  in^ose 
unfivorable  bias  on  the  outcome. 

Tayfim  (1980)  derived  the  probability  density  fimction  based  on  the  Stokes  profile 
approximation  up  to  the  second  order.  Huang  et  al.  (1983)  developed  the  probabihty 
density  fimction  considering  aU  terms  up  to  the  third  order.  These  probability  density 
fimctions  are  superior  to  the  Gram-Charlier  series  distribution  in  that  they  do  not  have  any 
negative  portion  over  the  entire  range  of  wave  displacement.  Con^arison  of  these 
probabihty  density  fimctions  and  wave  data  show  good  agreement. 

Langley  (1987)  apphed  a  probabihty  density  fimction  representing  the  response  of 
a  non-linear  system  developed  by  Kac  and  Siegert  (1947),  to  nonlinear  random  waves. 

This  probabihty  density  fimction  is  given  in  the  form  of  a  series  and  can  only  be  evaluated 
numericahy.  Ochi  and  Ahn  (1994a)  modified  the  Kac  Siegert  density  fimction  so  that  it  is 
presented  in  closed  form  with  parameters  determined  through  spectral  analysis.  They 
apphed  the  probabihty  density  fimction  to  non-Gaussian  coastal  waves.  In  the  apphcation 
process,  they  decomposed  the  spectra  into  linear  and  nonlinear  con:q)onents  to  clarify  the 
degree  of  non-linearity  involved  in  shahow  water  waves  (Ochi  and  Ahn  1994b).  The 
probabihty  density  fimction  does  not  have  any  negative  portion  in  its  entire  domain  and 
appears  to  represent  the  statistical  distribution  of  the  wave  profile  in  any  sea  state  during  a 
storm  in  deep,  intermediate  and  shahow  water. 
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Studies  concerned  with  the  non-Gaussian  properties  of  waves  in  finite  water  depth 
include  Thornton  and  Guza  (1983)  who  claimed  that  the  resuhs  of  their  analysis  of  field 
data  obtained  at  Torrey  Pines  Beach,  California,  show  that  the  Rayleigh  probability 
distribution  can  be  appUed  to  represent  wave  heights  observed  from  shallow  to  deep  (10.8 
meter)  water.  This  inches  that  all  observed  waves  can  be  considered  a  Gaussian  random 
process.  However,  the  sea  severity  in  relation  to  water  depth  at  all  locations  in  the 
nearshore  zone  where  their  measurements  were  taken  was  too  mild  and  below  the  limit  for 
waves  to  be  a  non-Gaussian  random  process.  This  subject  will  be  discussed  in  detail  in 
this  present  study. 

hi  closing  this  chapter,  h  is  inqiortant  to  note  that  there  still  remains  a  significant 
amount  of  work  to  do  in  refining  our  abihty  to  accurately  describe  irregular  seas  in  finite 
water  depths  and  over  a  wide  range  of  sea  severity. 


CHAPTERS 

BACKGROUND  /  STUDY  APPROACH 

3.1  Wave  Record  Data 

An  experiment  jointly  hosted  by  the  U.S.  Army  Coastal  Engineering  Research 
Center  (CERC)  and  the  National  Ocean  Survey  (NOS)  was  conducted  from  06  October 
1980  to  30  November  1980.  An  aspect  of  this  Atlantic  Ocean  Remote  Sensing  Land- 
Ocean  Experiment  (ARSLOE)  involved  ocean  wave  studies  in:q)lemented  near  the 
CERC’s  Field  Research  Facility  (FRF)  at  Duck,  North  Carolina.  In  situ  wave  measuring 
devices  were  dq)loyed  in  mean  water  depths  ranging  from  1.35  meters  to  24.4  meters  and 
were  placed  in  a  linear  array  extending  from  60  meters  offshore  to  12,000  meters  offshore. 
Figure  3.1. 

3.1.1  Data  Acquisition 

The  wave  displacement  data  analyzed  in  this  study  came  from  a  total  of  eight  wave 
measuring  devices.  The  FRF  has  a  550  meter  research  pier  oriented  perpendicular  to  the 
local  shore  which  was  equipped  with  seven  Baylor  resistance-type  wave  staffi.  These 
Baylor  wave  staffs  function  by  determining  the  change  in  electrical  resistance  experienced 
by  changes  in  water  elevations.  Round-trip  travel  time  of  an  electrical  current  from  an 
above- water  emitter  to  the  water  surface  is  measured.  Following  along  the  line  of  the 
pier  off-shore  was  a  Waverider  Buoy  Gauge.  This  is  a  surface-following  class  buoy.  It 
functions  by  double-integrating  the  output  of  an  accelerometer  held  vertical  by  a  pendulum 
to  produce  the  time  history  of  buoy  heave  motion.  The  statistics  of  the  buoy  heave  then 
represent  the  wave  displacements.  In  the  Waverider  buoy,  the  pendulmn  is  contained 
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entirely  within  the  spherical  float  so  the  buoy  does  not  have  to  be  stabilized.  These 
devices  recorded  both  wave  displacements  and  wave  periods  but  only  the  displacement 
data  were  analyzed  for  this  study.  The  wave  gauges  recorded  water  sur&ce  fluctuations 
simiiltaneously  at  all  eight  locations  at  a  frequency  of  4  Hertz  over  twenty  minute 
intervals.  This  resulted  in  4800  individual  wave  displacements  recorded  for  each  twenty 
minute  segments  continuously.  Due  to  the  significant  amount  of  data  obtained  during  this 
unique  experiment,  only  a  sub-set  of  the  data  was  analyzed.  This  chosen  data  sub-set 
began  at  1200  on  23  October  1980  and  ended  at  1315  on  25  October  1980.  This  two  day 

t 

segment  was  imique  in  that  a  severe  storm  occurred  from  23-25  October  •wiiich  was  the 
most  significant  meteorological  event  recorded  during  the  experiment.  Wind  q)eeds  on 
the  order  of  fifteen  meters  per  second  were  recorded  during  the  severe  storm  This  two 
day  sub-set  of  data  was  used  as  the  basis  for  this  study  since  it  included  the  entire  range  of 
wave  profiles  experienced  at  each  of  the  eight  wave  gauge  locations.  From  this  two  day 
sub-set,  a  total  of  nearly  5  million  individual  wave  displacements  were  analyzed  and  the 
wave  profile  data  covered  the  entire  range  of  sea  conditions  present  during  the  com|)lete 
56  day  ARSLOE  Project. 

3.1.2  Wave  Profile  Transformations 

The  data  collected  at  the  FRF  during  ARSLOE  is  unique  and  impressive.  It  is  very 
well  suited  for  analyzing  changes  in  the  wave  field  as  the  waves  propagate  through 
shoaling  depths.  The  bathymetry  at  the  experiment  location  is  relatively  siiiq)le  showing 
straight  and  parallel  contours.  This  provided  favorable  conditions  which  maximized 
control  over  the  environmental  impact  on  the  quality  of  the  data  received.  The  wave 
profile  data  provide  a  matrix  of  wave  profiles  allowing  both  spatial  and  teiiq)oral  analyses 
to  be  performed. 


11 


A  wide  range  of  wave  profiles  was  observed  during  the  e7q)erinient.  Figures  3.2- 
3.8  represent  a  mild  sea  condition  recorded  on  10/23  @  2155  for  all  finite  water  depth 
gage  locations  shown  in  Figure  3. 1.  Figures  3.9-3.15  represent  a  severe  sea  condition 
recorded  at  the  height  of  a  storm  on  10/25  @  1315  for  the  same  gage  locations.  It  is 
evident  fi'om  the  san:q)le  wave  data  that  the  wave  field  transitions  as  it  shoals  fi'om  deep  to 
shallow  water.  The  wave  data  representing  the  mild  sea  condition  at  Gage  625,  Figure 
3.8,  are  characterized  by  an  almost  symmetrical  distribution  of  positive  and  negative  wave 
profile  displacements.  This  is  characteristic  of  Gaussian  random  processes.  As  the  wave 
shoals  into  shallower  water,  the  wave  profile  transforms  showing  progressively  higher 
peaks  and  shallower  troughs.  This  process  is  clearly  non-Gaussian  and  must  be  analyzed 
differently.  It  is  also  clear  fi'om  observing  these  wave  records  that  the  degree  of  non¬ 
linearity  increases  as  the  relative  water  depth  decreases  and  as  sea  severity  increases.  This 
implies  that  analysis  of  this  nonlinearity  must  include  a  dependence  on  water  depth  as  it 
relates  to  the  magnitude  of  the  waves’  displacement  (sea  severity). 


Figure  3-3  Gage  635  - 10/23  @  2155  (d=2.83m) 


Figure  3-4  Gage  645  -  10/23  @  2155  (d=4.31in) 
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Figure  3-5  Gage  655  -  10/23  @  2155  (d=6.0m) 


Figure  3-6  Gage  665  -  10/23  @  2155  (d=6.08in) 


Figure  3-8  Gage  625  -  10/23  @  2155  (d=9.66m) 


Figure  3-9  Gage  615  - 10/25  @  1315  (d=2.23iii) 


Figure  3-11  Gage  645  -  10/25  @  1315  (d=4.66in) 


Figure  3-12  Gage  655  -  10/25  @  1315  (d=6.44m) 


Wyl/v  ^  v  wv  y  v  ‘y  V 


Figure  3-15  Gage  625  -  10/25  @  1315  (d=10.04m) 
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3.2  Nearshore  Zone  Characteristics 


All  wave  fields  propagating  in  finite  water  depths  are  influenced  by  the  attributes 
of  the  local  sea  environment.  These  attributes  include  water  depth  and  beach  slope.  An 
awareness  of  the  characteristics  of  these  attributes  is  in:q)ortant  to  understand  the  type  of 
influence  they  have  on  the  wave  field. 

3.2. 1  Discussion  of  Water  Deoth  Reeions  During  ARSLOE 


Based  on  linear  theory,  there  are  three  water  depth  regions  wfiich  are  defined  using 
the  velocity  potential  for  a  progressive  wave.  The  velocity  potential  is  defined  as 

4;rcosh(M)  T 

The  distinct  water  dq)th  regions  are  determiaed  by  the  asynq)totes  of  the  hyperbolic 
fimction.  Deep  water  wave  data  obtained  at  gage  location  710  and  shown  in  Table  3. 1 
allow  these  regions  to  be  defined  for  the  mild  sea  condition  and  the  severe  sea  condition. 

In  general,  the  deep  water  wave  condition  is  defined  as  kd  >  where  k  is  the  wave 
number  and  h  is  the  water  depth.  This  condition  occurs  at  d  =  21m  for  the  rmld  sea 
condition  on  10/23  @  2155  and  d  =  43m  for  the  severe  sea  condition  on  10/25  @  1315. 
The  shallow  water  condition  is  defined  as  kd  <  7t/10.  This  condition  occurs  at  d  =  0.6m  for 
the  mild  sea  condition  and  d  =  1.3m  for  the  severe  sea  condition  on  25  October.  The 
intermediate  water  condition  occurs  between  these  bormds.  All  wave  profile  data 
obtained  fi'om  ARSLOE  are  located  in  the  intermediate  water  depth  region. 

Table  3.1  Deep  Water  Wave  Statistics  for  MUd  and  Severe  Sea  Conditions 


Date/Time 

Zero  Crossing  Period,  T 

10/23  @  2155 

10/25  @  1315 

5.17  sec. 

7.43  sec. 

41.715  m 
86.157  m 
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3.2.2  Discussion  of  the  Surf  Zone  During  ARSLOE 

An  integral  aspect  of  any  study  wbich  covers  the  nearshore  region  nmst  analyze 
and  consider  the  surf  zone.  The  surf  zone  is  a  transient  entity  which  varies  in  width  off¬ 
shore  depending  on  the  local  sea  conditions.  These  local  conditions  include  the  beach 
profile,  sea  severity  and  direction  of  wave  field,  hi  order  to  define  the  surf  zone  in  this 
study,  wave  field  data  were  analyzed  vidiich  represented  the  range  of  sea  conditions 
experienced  over  the  two  day  period.  The  first  condition,  representing  the  mildest  sea 
condition,  was  taken  as  2155, 23  October  1980.  The  second  condition,  representing  the 
most  severe  sea  condition  was  taken  as  1315,  25  October  1980. 

Linear  wave  theory  will  be  used  to  determine  the  effects  of  shoaling  as  the  wave  field 
enters  shallow  water.  Refi'action  will  not  be  considered  since  this  process  tends  to  reduce 
the  shoaling  wave  heights  and  thus  move  the  surf  zone  shoreward.  This  is  shown  by  the 
definition  of  the  refi’action  coefficient,  where  0o  and  0  are  the  deep  water 

V  cos^ 

wave  angle  and  the  finite  water  depth  wave  angle  respectively.  This  ratio  will  always  be 
less  than  one  causing  a  reduction  in  the  resulting  wave  height.  In  order  to  ensure  that 
conservative  results  are  obtained,  refraction  will  be  neglected  in  this  study.  The  beach 
profile  at  the  ARSLOE  project  site  shows  slightly  varying  beach  slopes  so  a  wave 
breaking  criteria  which  considers  the  beach  slope  will  be  used.  Weggel  (1972) 
reinterpreted  laboratory  results  and  showed  a  dependency  of  breaker  height  on  beach 
slope.  The  following  equations  were  used  to  define  the  surf  zone  locations; 


= 


g  1 


g2 


Conservation  of  Energy;  H  ^ 
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(assumed  constant),  g  =  gravitational  acceleration,  k  =  the  wave  number  and  d  =  local 
water  depth. 

hi  order  to  standardize  the  location  of  the  surf  zone,  the  root-mean-square  wave 
height,  H„ns ,  will  be  used  to  represent  the  irregular  wave  field.  H  = 

vs^ere  Hi  is  defined  as  the  wave  height  of  the  i***  wave  of  the  irregular  wave  field.  The  surf 
zone  which  results  from  the  use  of  Hms  will  represent  the  approximate  “center  of  mass” 
location  where  the  most  significant  set  of  individual  waves  in  the  wave  field  begin  to 
break.  The  surf  zone  calculations  were  based  on  a  root-mean-square  wave  height  that  was 
determined  from  the  deep  water  wave  displacement  data.  This  single  wave  height  was 
used  to  represent  the  irregular  sea  and  shoaled  landward.  Since  the  real  sea  will  have  a 
significant  number  of  waves  with  heights  larger  and  smaller  than  the  root-mean-square,  the 
actual  breaking  zone  is  an  obscure  entity.  This  technique  for  determining  the  surf  zone 
locations  is  adequate  to  make  general  conqiarisons  when  analyzing  the  spatial  changes  in 
the  nonlinearity  of  the  irregular  sea  in  finite  water  depths.  The  surf  zone  locations  during 
the  mildest  and  most  severe  sea  conditions  over  the  two  day  period  are  shown  in  Table 
3.2. 

In  Table  3.2  it  is  shown  that  during  the  mild  sea  condition,  the  breakpoint,  indicating 
the  beginning  of  the  surf  zone,  is  located  70  meters  from  the  shoreline  in  1.5 1  meters  of 
water.  This  brealq)oint  is  the  location  where  a  wave  having  a  height  equivalent  to  the 
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root-mean-square  wave  height  of  the  wave  field  will  break.  During  the  severe  sea 
condition,  this  breakpoint  occurs  110  meters  from  the  shoreline  in  3.32  meters  of  water. 
Realistically,  waves  would  be  observed  breaking  considerably  seaward  of  the  breakpoints 
determined  in  the  respective  sea  conditions  due  to  the  spectrum  of  wave  heights  present  in 
an  irregular  sea.  The  root-mean-square  wave  height  breakpoint  is  used  in  order  to 
standardize  the  position  of  the  surf  zone  for  canq)arison  purposes. 


Table  3.2  Surf  Zone  Locations  for  Mild  and  Severe  Sea  Conditions 


Sea  Condition 

Start  of  Surf  Zone 

height  Hniisy  Breaks 

-Dist.  from  Shore- 

1.51  m 

70  m 

3.32  m 

110m 

3.3  Non-Gaussian  Probabihtv  Density  Function 

The  probabilistic  approach  used  in  this  study  will  make  use  of  the  probability 
density  function  apphcable  to  non-gaussian  random  processes  developed  by  Ochi  and  Ahn 
(1994a).  This  probability  density  function  appears  to  represent  non-Gaussian  waves  in 
shallow  and  intermediate  water  depths  as  well  as  Gaussian  waves  in  deep  water. 

3.3.1  Non-Gaussian  Probability  Density  Function  Background 

The  probability  density  function  has  a  background  of  the  Kac-Siegert  solution 
developed  for  the  response  of  non-linear  systems.  The  density  function,  y(t),  is  given  as 
the  sum  of  the  standardized  normal  variates  and  its  squares  as  follows 
y(t)=  + 
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The  parameters  Pj  and  are  determined  by  finding  the  eigenfimction  and  eigenvalues 
of  the  integral  equation  given  by 

where  K((Di  ,02)  =  H(oi ,  CO2  )  yfS(^^)S(&^ 

\|/j(©)  =  orthogonal  eigenfimction  =  (  J j(a))y/ ^(k)do)  =  1  for  j  =  k 

otherwise,  v|/j(©)  =  0 

H(©i  ,(02)  =  second  order  fi'equency  response  fimction 
S(©)  =  output  spectral  density  fimction 

As  can  be  seen  fi-om  the  form  of  this  solution,  the  probability  density  fimction  of  y(t), 
can  not  be  found  without  knowledge  of  the  second  order  fi’equency  response  fimction, 
H(©i ,  ©2 )  and  there  is  no  way  to  evaluate  this  response  fimction  for  irregular  seas.  Ochi 
and  Ahn  showed  however,  that  Kac-Siegerts’  solution  can  be  obtained  fiom  the  time 
history  of  the  wave  record  without  any  knowledge  of  the  second  order  fiequency  response 
fimction,  H(©i ,  ©2 ).  The  probability  density  fimction  thus  obtained  for  describing  non- 
gaussian  random  processes  is 

=  -7=i — exp  -  — — — — (1  -  yoM.  -  exp(-rax)y  -  yax 
•sjlna,  L  lyyaa^Y 

,  L28^r;x>0 
where,  /  =  1 

^  3.00/or:x<0 

The  parameters,  a,  |i»  and  are  determined  fiom  the  cumulants  evaluated  fiom  the  wave 
data.  These  cumulants  are  related  to  the  Pj  and  Aj  values  of  the  Kac-Siegert  solution  by 
the  following  relationships 

2  A' 

k,  =  s  . 

J  =  1 

2  IN 

ks  =  6 ^  jX  .  +  8 ^  X^ j  ,  etc. 

7  =  1  3=\ 
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In  the  present  study,  ki  represents  the  mean  water  level  and  is  therefore  defined  by, 
ki  =  0.  k2  and  ks  are  equal  to  the  data  sanq)le  moments,  E[x^]  and  E[x^]  respectivefy. 
Based  on  this  relationship,  the  parameters,  a,  p*  and  can  be  determined  through  the 
iterative  solution  of  three  equations  along  with  the  analysis  of  a  data  record  for  the  non- 
Gaussian  random  process  being  analyzed.  The  three  equations  are 
a  +  a  +  p,  =  E[x]  =  0 
a^-2aV.  =  E[x^] 

2aaV3-  8aV*)  =  E[x^] 


3.3.2  Probabilitv  Density  Function  Application  to  Wave  Data 

In  order  to  determine  values  for  the  three  parameters,  a,  p»  and  ,  which  are 
used  to  define  the  degree  of  non-linearity  a  shoaling  wave  field  has  at  a  given  location  and 
sea  state,  data  records  must  be  analyzed.  With  the  applicabihty  of  the  probability  density 
function  verified  based  on  the  wave  data  collected  during  ARSLOE,  the  approximately 
1000  data  records  fi’om  1200  on  23  October  to  1315  on  25  October  will  be  analyzed  to 
determine  values  for  the  three  parameters  defining  the  probabihty  density  function.  The 
procedure  used  during  this  study  to  perform  this  task  is  as  follows: 

-  Compute  the  standard  deviation  of  each  twenty  minute  wave  record  consisting  of 

4800  discrete  wave  displacement  data  points  E{x^ )  =  “  T!  • 

(E[x]  does  not  need  to  be  computed  since  it  is  defined  as  zero  and  the  data  are 

aheady  standardized). 

-  Compute  the  third  moment  of  the  same  twenty  minute  wave  record 

rii 
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-  Define  the  standard  deviation  of  the  data  record  as  cr  =  ylE(x^)  . 

-  Evaluate  the  three  parameters,  a,  and  ,  through  iterative  solution  of  the 

following  three  equations; 

a  +  a  +  p.  =  E[x]  =  0 
a^-2aV.  =  E[x^] 

2acjV3-8aV.)  =  E[x^]. 

It  is  important  to  note  that  the  standard  deviation  of  the  data  record  is  used  in  the 
analysis  as  an  indicator  of  the  sea  severity  and  not  the  significant  wave  height.  This  is 
because  the  definition  of  significant  wave  height  is  not  applicable  to  non-Gaussian  waves. 

All  available  data  sets  covering  the  seven  finite  water  depth  gage  locations 
between  1200,  23  October  1980  and  1400,  25  October  1980  were  analyzed  and  a,  p»  and 
were  determined. 

Based  on  a  plot  of  the  individual  parameters  in  the  space  domain,  it  was 
determined  that  a  smoothing  technique  should  be  used  to  refine  the  “a”  parameter  values. 
This  smoothing  was  done  in  order  to  remove  the  sporadic  tertiary  effects  while  ensuring 
the  smoothed  “a”  value  does  not  lose  any  authenticity.  An  e7q)onential  smoothing 
techmque  was  used  which  resulted  in  the  lowest  mean  absolute  deviation  (MAD)  between 
the  raw  “a”  parameter  value  and  the  smoothed  “a”  parameter.  The  remaining  two 
parameters,  p»  and  did  not  show  the  same  susceptibility  to  tertiary  effects  so  no 

smoothing  was  done.  Comparisons  between  the  actual  and  smoothed  “a”  parameters  are 
shown  in  Appendix  A  for  all  Gage  positions. 

3.3.3  Dimensional  Analysis  of  PDF  Parameters 

For  each  of  the  approximately  one-hundred-fifty  data  sets  per  gage  position,  the 
probability  density  fimction  parameters,  a,  p. ,  a*  were  determined  iteratively.  Each  of 
the  parameters  contain  discrete  dimensions.  Based  on  dimen.ginnal  analysis,  these 
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dimensions  are  as  follows:  It  is  known  that  k2,  representing  the  data  sets’  variance  has 
dimensions  of  [L^].  It  is  also  known  that  ks,  representing  the  data  sets’  third  moment  has 
dimensions  of  [L^].  It  follows  from  the  knowledge  of  these  dimensions  that  the 
parameters,  a,  p* ,  c*  must  have  dimensions  [1/L],  [L]  and  [L]  req)ectively. 

In  order  to  minimize  any  dependence  on  external  conditions  and  to  ensure  the 
solutions  are  universally  applicable,  these  parameters,  a,  p» ,  a*  must  be  made  non- 
dimensional  The  non-dimensional  quantities  that  Avill  be  used  in  this  study  are,  aa,  p*/cT, 
aja  and  a/d.  a  is  defined  as  the  standard  deviation  of  that  wave  displacement  data  and  d 
is  defined  as  the  water  depth. 


CHAPTER  4 

PROBABILnY  DENSITY  FUNCTION  ANALYSIS 
4.1  Probability  Densitv  Function  Applicabilitv  Verification 

As  stated  in  the  introduction,  one  of  the  objectives  of  this  study  is  to  determine 
vdiether  or  not  the  probability  density  Sanction  developed  by  Ochi  and  Ahn  is  valid  for 
representing  waves  in  finite  water  depths.  To  verify  the  applicability  of  this  probability 
density  fimction,  wave  displacement  data  sets  obtained  fi^om  the  ARSLOE  eaqperiment 
were  used  in  conqaaring  the  histogram  of  the  data  record  to  the  probability  density 
fimction.  Each  data  set  contained  4800  wave  displacement  readings  taken  over  a  twenty 
minute  period.  To  ensure  a  thorough  verification  of  the  probability  density  fimction,  the 
data  records  used  for  the  analysis  were  fi'om  all  seven  “finite”  water  depth  gage  locations 
and  analysis  was  done  over  a  wide  range  of  sea  state  conditions.  One  analysis  was  done 
on  data  firom  23  October  1980  at  2155.  This  represented  a  mild  sea  state.  The  other 
analysis  was  done  on  data  fi'om  25  October  1980  at  1315.  This  data  set  was  at  the  peak 
of  the  most  severe  storm  ejqjerienced  during  the  entire  ARSLOE  study. 

Figures  4. 1-4. 14  show  comparisons  between  the  histogram  data  and  the  non- 
Gaussian  probability  density  fimction.  Included  also  in  the  figures  is  the  Gaussian 
(normal)  probabihty  density  fimction  for  comparison.  Figures  4. 1-4.7  pertain  to  the  mild 
sea  condition  \\Me  Figures  4.8-4. 14  represent  the  severe  sea  condition.  As  can  be  seen 
fiom  these  figures,  the  probability  density  fimction  agrees  very  well  with  the  histogram 
generated  fiom  the  data  record  in  all  cases.  This  imphes  that  the  probability  density 
fimction  derived  by  Ochi  and  Ahn  can  be  used  with  confidence  when  analyzing  non- 
Gaussian  random  processes  such  as  wave  displacements  in  shallow  and  intermediate 
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waters.  It  can  also  be  seen  that  the  probability  density  function  approximates  the  normal 
distribution  wiien  the  wave  field  is  in  the  deeper  intermediate  water  condition  under  mild 
sea  conditions. 


Gage  615-10/23  @  2155 


Displac«mtnt  (tn) 


His  to  gram  Data  - PDF . Gaussian 


Table  4.1  Gage  615  -  10/23  @  2155 


WAVE  PROPERTIES 
615  10/23(5)2155 

a  (l/m)  .572645 


luCni) 

CT.(ni) 

CT(in) 

d(m) 

a/d 


- .04948 
.275683 
.267219 
1.74 
.153574 


Figure  4.1  Histogram  Data 


Probability  Density  Function  Probability  Density 
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Table  4.2  Gage  635  -  10/23  @  2155 


WAVE  PROPERTIES 


10/23(2)2155 


.481441 


- .0547 


.318871 


.30989 


.109502 


Table  4.3  Gage  645  -  10/23  @  2155 


WAVE  PROPERTIES 


23(2)2155 


.200988 


.03084 


.343446 


.340708 


.079051 


Figure  4.3  Histogram  Data 
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Oaae  666-10/23^2166 


1.2 
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Figure  4.4  Histogram  Data 


Gage  666  - 10/23  @  2166 


1.2 
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Table  4.4  Gage  655  - 10/23  @  2155 


WAVE  PROPERTIES 

655 

10/23(2),  2155 

a  (l/m) 

.104437 

fi.  (m) 

- .01771 

<T.(lll) 

.340056 

tr(m) 

.339137 

d(ni) 

6.0 

CT/d 

.056523 

Table  4.5  Gage  665  -  10/23  @  2155 


WAVE  PROPERTIES 

665 

10/23  (2).  2155 

a  (1/m) 

.126266 

P.  (m) 

-.01835 

<T.(m) 

.365155 

a(m) 

.364236 

d(m) 

6.08 

cr/d 

.059907 

Figure  4.5  Histogram  Data 
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Gage  ^6 -10/23^2166 
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Table  4.6  Gage  675  - 10/23  @  2155 


WAVE  PROPERTIES 

67S 

10/23  (2),  2155 

a  (l/m) 

.094349 

(m) 

- .01014 

0.(01) 

.375221 

o(di) 

,374948 

d(ni) 

6.82 

o/d 

,054978 

Figure  4.6  Histogram  Data 


Gage  626 -10/23  @2166 


Displacement  (m) 


EE2E2223  Hstoyam  Data  - PDF  . Gaussian 


Table  4.7  Gage  625  -  10/23  @  2155 


WAVE  PROPERTIES 

625 

a  (1/m) 

.069914 

P.  (m) 

-.00819 

<T,(m) 

.394571 

<T(m) 

.394401 

d(m) 

9.66 

a/d 

.040828 

Figure  4.7  Histogram  Data 
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Gage  615-10/25  ®  1315 


Disphie»m*iit  (m) 


vjuwMitJi  Histogram  Data  - PDF . Gaussian 


Table  4.8  Gage  615  - 10/25  @  1315 


WAVE  PROPERTIES 

615 

10/25  (®,  1315 

a  (1/m) 

.372027 

H.(m) 

-.09863 

o,(m) 

.498230 

CT(m) 

.47981 

d(m) 

2.23 

<7/d 

.215161 

Figure  4.8  Histogram  Data 


Gage  636  - 10/26  @  1316 


Displacement  (m) 


Pflgggggj  Hstogram  Data  - PDF  . Gaussian 


Table  4.9  Gage  635  - 10/25  @  1315 


WAVE  PROPERTIES 

635 

10/25  (Si  1315 

a  (1/m) 

.318011 

P.  (m) 

- .10497 

a.  (m) 

.591069 

CT(m) 

.573397 

d(m) 

3.32 

a/d 

.17271 

Figure  4.9  Histogram  Data 


Probability  Density  Function 
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Figure  4.10  Histogram  Data 


Table  4.11  Gage  655  -10/25  @  1315 
I  WAVE  PROPERTIES 


a  (l/m) 

.269927 

- .12356 

a.(in) 

.658813 

CT(m) 

.636821 

d(m) 

6.44 

a/d 

.098885 

Figure  4.11  Histogram  Data 
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Gage  665  -  1(tf25  e  1315 


0.8 


Displacement  (m) 


^^^^Histo^'amData  - PDF  . GaiBsian 


Table  4.12  Gage  665  -10/25  @  1315 


WAVE  PROPERTIES 

66S 

10/25  (S),  1315 

a  (l/m) 

.269743 

- .12978 

CT.(m) 

.67434 

CT(m) 

.650698 

ddn) 

6.52 

o/d 

.09980 

Figure  4.12  Histogram  Data 


Gage  676  - 10/26  e  1316 


Table  4.13  Gage  675  -10/25  @  1315 


WAVE  PROPERTIES 

675 

10/25  ®  1315 

a  (l/m) 

.220143 

-.12062 

CT,(m) 

.727669 

<T(m) 

.708477 

d(m) 

7.22 

a/d 

.098127 

Figure  4.13  Histogram  Data 
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Table  4.14  Gage  625  -10/25  @  1315 

WAVE  PROPERTIES 

625  10/25  (a).  1315 

a  (l/m) 

.140009 

H.  (m) 

- .10961 

CT.  (m) 

.872244 

®  (m) 

.858792 

d(iii) 

10.04 

cr/d 

.085537 

Figure  4.14  EKstogram  Data 


Based  on  a  comparison  with  the  water  depth  regions  defined  in  Section  3.2. 1,  it 
can  be  seen  that  the  histogram  data  and  the  probahihty  density  fimction  closely  resemble 
the  normal  distribution  even  after  the  wave  field  enters  intermediate  water  for  milder  sea 
conditions.  This  is  shown  in  Figures  4.6  and  4.7.  This  in:q)lies  that  the  characteristics  of 
the  wave  profile  are  dependent  on  the  sea  condition  as  well  as  the  water  depth. 

The  histogram  data  for  both  the  mild  and  severe  sea  conditions  in  Figures  4. 1-4. 14 
show  the  distribution  of  the  wave  profiles  progressively  becomes  more  non-Gaussian  as 
the  waves  propagate  toward  shallower  depths.  The  non-Gaussian  trend  becomes  more 
extreme  as  the  sea  condition  becomes  more  severe.  It  is  significant  to  note  also  that  the 
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non-Gaussian  characteristics  are  most  pronounced  at  the  shallower  gage  locations,  Gages 
615  and  635. 


4.2  Gaussian/non-Gaussian  Boundary  Definition 

The  histogram  data  shown  in  Figures  4.1-4. 14  represent  the  same  wave  profiles 

depicted  in  Figures  3.2  -3.15.  The  transformation  from  Gaussian  to  non-Gaussian 

profiles  begin  with  subtle  changes  in  the  anq)litude  of  the  negative  displacements  and  the 

width  of  the  troughs.  This  initial  transformation  can  be  seen  be  con^aring  Figures  4.4  and 

» 

4.6.  Figure  4.6  shows  a  nearly  normal  distribution.  Figure  4.4  ^ows  a  slight  elevation  of 
the  apex  and  a  narrowing  of  the  negative  tail.  It  is  also  significant  to  note  that  the  positive 
part  of  the  distribution  remains  relatively  unchanged.  This  indicates  that  the  initial  changes 
affect  the  trough  region  of  the  profile  more  drastically  than  the  positive  region.  Figure  4.4 
demonstrates  that  the  trough  region  becomes  ^allower  and  wider.  As  the  profile 
continues  to  ^oal,  changes  begin  to  occm  in  the  positive  displacement  region  of  the 
distribution.  The  peaks  become  narrower  and  higher.  This  can  be  seen  graphically  by 
con^aring  Figures  4.2  and  4.4.  The  positive  region  of  the  wave  profile  is  not  as  full  when 
cony)ared  with  the  normal  distribution  and  the  positive  tail  extends  farther  than  the  normal 
curve. 

As  described  above,  significant  transformations  occur  in  the  wave  profile  as  the 
wave  train  propagates  into  shallower  water.  There  is  a  point  in  this  transformation  when 
the  wave  profile  deviates  significantly  enough  from  the  normal  distribution  that  the  waves 
can  no  longer  be  considered  a  Gaussian  random  process.  Identification  of  this  boundary  is 
mq)ortant  to  ensure  assumptions  made  regarding  the  waves  remain  valid.  Analysis  of  the 
histogram  data  shows  that  the  non-Gaussian  characteristics  of  coastal  waves  depends  on 
water  depth  and  sea  severity.  Based  initially  on  analysis  of  the  histogram  data  shown  in 
Figures  4.5  and  4.6,  it  appears  this  transition  boundary  occurs  when  the  non-dimensional. 
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wave  standard  deviation,  o/d,  is  approximately  0.06.  The  non-dimensional  parameter, 
a/d,  is  used  because  it  considers  both  sea  severity  and  water  depth.  Eight  random  sea 
conditions  were  selected  having  a/d  values  between  0.04  and  0. 10.  These  eight  data  sets 
are  depicted  in  Figure  4. 15  in  terms  of  aa  as  a  function  of  a/d .  The  data  points  labeled  A- 
H  in  Figure  4.23,  each  correspond  to  a  distribution  shown  in  Figures  4. 15-4.23.  Data 
points  A-G  were  intentionally  grouped  around  the  10  meter  water  depth.  This  minimized 
any  inqiact  water  depth  and/or  beach  slope  may  have  on  the  analysis  results.  The 
statistical  data  for  these  points  are  provided  in  Table  4. 15. 
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Figure  4.15  Sample  Data  Points  to  Analyze  Gaussian/Non-Gaussian  Boundary 
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Table  4.15  Statistical  Properties  of  Data  Points  A-H  of  Figure  4.15 


Figure  4.16  Point  A  of  Figure  4.15 
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Comparison  with  Norma i  Distribution 
aid  =  .080 


Prob.  Density  Fen. . Normal  Distribution 


F^ure  4.17  Point  B  of  Figure  4.15 


Comparison  with  Norma i  Distribution 
old  =  .071 


Rob.  Density  Fen. . Normal  Distribution 


Figure  4.18  Point  C  of  Figure  4.15 


37 


Figure  4.19  Point  D  of  Figure  4.15 
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Comparison  with  Normai  Distribution 
aid  -  .062 


Rrob.  Density  Fen. . Normal  Dstribution 


Figure  4.21  Point  F  of  F^ure  4.15 


Comparison  with  Normal  Distribution 
aid  =  .0398 


FYob.  Density  Fen. . Normal  Distributton 


Figure  4.22  Point  G  of  F^ure  4.15 
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Comparison  with  Normal  Distribution 
a/d  = .069 


Displacement  (m) 

- FYob.  Density  Fen. . Nornral  Distributton 

Figure  4.23  Point  H  of  Figure  4.15 

Figures  4.16  -  4.23  confirm  that  the  region  between  a/d  =  0.06  and  a/d  =  0.08, 
shown  in  Figures  4.20  and  4. 17  respectively,  appears  to  be  when  the  wave  profile 
transforms  fi’om  Gaussian  to  non-Gaussian.  A  closer  inspection  of  this  region  is  shown  in 
Figures  4. 18  -  4.20.  The  distribution  shows  a  considerable  increase  in  deviation  fi'om  the 
normal  between  a/d  =  0.06  and  a/d  =  0.0616  as  shown  in  Figures  4.20  and  4. 19 
respectively. 

Due  to  the  subjective  nature  of  this  boimdary,  a  general  boundary  can  be  defined 
based  on  probabihstic  analysis.  This  boundary,  defining  the  conditions  under  which  the 
wave  profile  can  no  longer  be  considered  Gaussian,  occurs  when  a/d  =0.06. 

A  similar  study  was  conducted  by  Ochi  and  Wang  (1984);  however,  they  defined 
the  boundary  condition  between  Gaussian  and  non-Gaussian  wave  profiles  in  terms  of  the 
wave  record  parameter,  X3  represents  the  skewness  of  the  wave  profile  which  they 

determined  to  be  the  dominant  parameter  affecting  the  non-Gaussian  characteristics  of 
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coastal  waves.  Ochi  and  Wang  concluded  that  coastal  waves  in  seas  for  wdiich  X3  is  less 
than  0.2  can  be  considered  as  a  Gaussian  random  process.  For  con:q)arison  purposes,  this 
boundary  condition  is  plotted  along  with  the  findings  of  this  present  study  in  Figure  4.24. 
In  terms  con^atible  with  the  approach  of  this  present  study,  their  boundary  condition  can 
be  defined  as  a  /d  =  0.0543.  It  is  encouraging  that  there  shows  very  close  agreement 
between  both  findings.  Additional  analysis  was  contributed  by  Ochi  and  Ahn  (1994b) 
whose  data  points  are  shown  in  Figure  4.24  as  white  circles.  Ochi  and  Ahn  (1994b) 
defined  the  transition  criterion  between  Gaussian  and  non-Gaussian  as  occurring  \^^en  a  < 
0.01  w^ere  “a”  is  one  of  the  parameters  in  the  non-Gaussian  probability  density  fimction. 
The  white  circle  data  points  define  the  limiting  cases  based  on  his  criterion.  This  third 
method  for  determining  the  boundary  between  Gaussian  and  non-Gaussian  wave  profiles 
agrees  very  well  with  the  other  two  fiirther  veriJfying  the  validity  of  this  finding. 
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Figure  4.24  Gaussian/Non-Gaussian  Boundary  Deflned 


CHAPTERS 

TREND  ANALYSIS  OF  PARAMETERS 


5.1  Relationship  between  a  and  Water  Depth 

As  was  discussed  earlier,  the  distribution  of  the  wave  profile  is  highly  dependent 
on  sea  severity  and  water  depth  with  sea  severity  being  defined  in  terms  of  the  wave 
profiles’  standard  deviation,  a.'  Since  these  two  terms  will  be  used  considerably  in  the 
following  trend  analyses,  their  relationsh^  with  each  other  will  be  examined. 

Figure  5. 1  shows  examples  of  simultaneous  measurements  of  ct  plotted  as  a 
function  of  distance  offihore  for  two  cases.  It  is  evident  fi'om  Figure  5. 1,  that  the 
standard  deviation  of  the  wave  profile  is  influenced  by  water  depth.  As  water  depth 
decreases,  the  relative  magnitude  of  the  standard  deviation  decreases.  It  is  also  seen  that 
for  more  severe  sea  conditions,  represented  by  10/25  @  1315  data  in  Figure  5.1,  the 
overall  change  in  0  as  the  wave  field  propagates  into  shallower  waters  is  greater  as 
con^)ared  with  milder  sea  conditions  (10/23  @  2155).  Based  on  theoretical 
considerations  of  0,  it  can  be  assumed  that  0  approaches  zero  as  distance  offshore 
approaches  zero.  This  trend  is  observed  in  Figure  5. 1  since  both  cases  shown  appear  to 
converge  to  the  origin.  From  the  information  on  the  surf  zone  locations  for  these  two 
cases  confuted  in  Chapter  3,  it  appears  that  0  begins  to  decrease  rapidly  towards  zero 
once  the  break  point  is  reached. 

The  information  shown  in  Figure  5. 1  based  on  two  cases  is  shown  in  Figure  5.2  for 
all  data  sets  analyzed  in  this  study.  Instead  of  0  being  plotted  as  a  function  of  distance 
offshore.  Figure  5.2  plots  0  as  a  function  of  water  depth.  This  will  take  into  consideration 
any  local  changes  in  mean  sea  level. 
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Figure  5.1  Spatial  Comparison  of  o  at  Mild  and  Severe  Sea  Conditions 


The  trends  evident  in  Figure  5.2  reinforce  the  analysis  done  for  the  two  cases  in 
Figure  5. 1.  The  upper  bounds  of  the  discrete  data  rq)resent  the  severe  sea  conditions 
while  the  lower  bounds  represents  the  mUd  sea  conditions. 
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It  can  also  be  inferred  from  Figure  5. 1  and  Figure  5.2  that  as  the  sea  state  becomes 
milder,  the  waves’  standard  deviation  approaches  a  constant  value  and  hence  the  sea 
severity  approaches  a  constant  before  decreasing  towards  the  origin. 

hi  the  absence  of  elaborate  wave  data,  the  value  for  o  can  be  readily  determined 
with  a  knowledge  of  the  local  and  breaking  wave  heights  and  the  relative  water  depth, 
defined  by  — ^  ,  where  g  =  gravitational  constant  and  T  =  wave  period,  using  a  method 

such  as  stream  fimction  theory  (Dean  1974).  To  ensure  the  use  of  accurate  wave  height 
values  in  applying  the  stream  fimction  theory,  nonlinear  shoaling  methods  should  be  used 
in  finite  water  depths  to  determine  the  breaking  and  local  wave  heists. 

5.2  Non-Gaussian  PDF  Parameter  Analysis 

The  non-Gaussian  probability  density  fimction  parameters,  a,  p*  and  ct»,  will  be 
analyzed  in  non-dimensional  form  spatially  as  a  fimction  of  distance  ofishore,  against  each 
other  in  non-dimensional  form,  and  with  aa  as  a  fimction  of  a/d.  The  first  analysis  of  the 
non-Gaussian  probability  density  fimction  parameters  will  involve  coixqparing  the  non- 
dimensional  parameters  against  the  distance  the  gage  was  located  offrhore.  This  will 
provide  a  spatial  orientation  of  the  parameters  which  will  disregard  any  local  fluctuations 
in  water  depth  due  to  tidal  cycle  of  storm  surge.  The  results  are  shown  in  Figures  5. 3-5. 5. 

5.2.1  Spatial  Analysis  of  Parameters 

All  three  parameters  shown  in  Figures  5. 3- 5. 5  appear  to  gradually  increase  in 
absolute  magnitude  until  approximately  150  meters  to  175  meters  offshore.  The  mild 
condition  represented  by  10/23  @2155  consistently  yields  lower  values  when  cotr^ared 
with  the  more  severe  sea  condition  on  10/25  @1315.  As  the  sea  condition  becomes  more 
mild,  the  wave  profile  resembles  the  normal  distribution.  When  the  wave  profile 
represents  exactly  the  normal  distribution,  “a”  equals  zero  and  o*  equals  a.  It  can  be  seen 


Plot  of  Bottom  Profile  and  "aa"  Comparison 
Distance  Offshore  (m) 
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Figure  5.3  Spatial  Comparison  of  ao  at  Mild  and  Severe  Sea  Conditions 


Plot  of  Bottom  Profile  and  "a*  /a"  Comparison 
Distance  Offshore  (m) 


Figure  5.5  Spatial  Comparison  of  o./ct  at  Mild  and  Severe  Sea  Conditions 
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from  Figures  5.3  and  5.5  that  as  approaches  zero  under  mild  conditions  and  s*/s 
approaches  one.  Figure  5.4  shows  that  m*/s  approaches  zero  under  these  same  conditions 
\^ch  is  also  consistent  with  the  theoretical  limit  of  the  non-Gaussian  probability  density 
fimction. 

There  is  also  a  significant  perturbation  evident  in  Figures  5. 3-5. 5  around  150 
meters  offihore.  The  mild  sea  condition  shows  this  abrupt  change  occurring  seaward 
v^en  con^ared  with  its  location  under  severe  sea  conditions  and  the  perturbation  is  more 
pronoimced  in  the  severe  sea  condition  case.  In  general,  all  trends  analyzed  for  the  non- 
dimensional  parameter,  m*/s,  are  almost  identical  results  to  those  fomd  for  as  and  s*/s  but 
in  the  mirror  image.  This  is  due  to  the  fact  that  the  value  of  m>  is  negative.  The  negative 
sign  is  supported  by  the  theoiy  used  to  develop  the  probability  density  functions’  three 
governing  equations. 

The  next  spatial  analysis  will  consist  of  conq)aring  the  non-dimensional  parameters 
with  water  depth.  This  conq)arison  will  be  sensitive  to  any  influence  of  local  depth 
changes  occurring  from  tidal  cycles  and/or  storm  surges.  Figures  5.6  -  5.8  show  the 
results  of  this  analysis. 

The  trends  from  this  analysis  are  the  same  as  those  determined  by  Figures  5.3-5. 5. 
Due  to  the  difference  in  scale  of  the  abscissa  the  curves  based  on  water  depth  appear  to 
show  much  smoother  transitions.  Figures  5. 6-5. 8  include  two  additional  cases 
representing  intermediate  sea  conditions  between  those  represented  by  10/23  @  2155  and 
10/25  @  1315.  These  two  additional  cases  approximate  the  trends  displayed  by  the  mild 
sea  condition  of  10/23  @2155.  This  is  consistent  with  the  data  and  time  of  these  cases 
since  they  occiirred  prior  to  the  severe  storm  so  the  sea  condition  was  not  yet  inq)acted 
significantly. 

As  was  discussed  earher,  in  the  deeper  water  conditions  all  three  parameters  tend 
to  show  convergence  towards  their  respective  values  for  a  Gaussian  wave  profile.  In 
extremely  shallow  water  conditions,  the  parameters  also  tend  to  converge  towards  a 
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Figure  5.8  o./ct  as  a  Function  of  Water  Depth 
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constant  value  wiiich  is  different  depending  on  the  sea  condition.  For  the  milder  sea 
conditions,  the  parameters  tend  to  converge  to  the  values  shown  in  Table  5.1. 


Table  5.1  General  Shallow  Water  Limits  of  Parameters 


Non-Dimensional  Parameter 

Mild  -  10/23  ®  2155 

Severe  -  10/25  (%  1315 

aa 

0.15 

0.18 

li./a 

-0.18 

-0.22  to  -0.25 

a*/a 

1.033 

1.04  to  1.045 

When  the  separate  plots  for  each  non-dimensional  parameter  are  con^ared  in  both 
spatial  cases,  it  is  apparent  that  similar  trend  descriptions  fit  all  three  parameters.  For  this 
reason,  a  closer  analysis  will  be  made  among  the  three  variables  to  determine  if  this 
similarity  can  be  quantified.  The  next  trend  analysis  consists  of  con^aring  the  three 
parameters  with  each  other. 

5.2.2  Relationshy  Between  the  three  non-dimensional  parameters,  ag.  and  gJo 

As  discussed  earlier,  the  three  non-dimensional  parameters  appear  to  be  highly 
correlated  with  each  other.  A  closer  look  will  be  taken  to  see  if  specific  relations  can  be 
developed  connecting  the  three  non-dimensional  parameters  together.  Figure  5.9  and 
Figure  5. 10  show  plots  of  the  raw  data  representing  \x,Jo  and  a*/a  conqiared  with  aa 
respectively  at  the  extreme  water  depth  locations,  gage  615  and  gage  625. 

All  data  sets  jfrom  Gages  615  and  625  were  plotted  since  they  represent  the 
shallowest  and  deepest  gage  locations.  There  are  very  definite  relationships  which  hold 
constant  between  the  extreme  range  of  water  depths  and  sea  conditions,  that  link  ao  to 
|a.»/o  and  a»/a.  The  fimctional  relationship  can  be  determined  only  after  considering  all 
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the  data  points  from  all  water  depth  and  sea  severity  conditions  to  ensure  it  remains 
consistently  valid.  A  deterministic  relationship  between  these  non-dimensional  parameters 
is  a  topic  of  active  research. 


Figure  5.9  Plot  of  p,»/cj  as  a  function  of  aa  for  Gages  615  and  625 


Figure  5. 10  Plot  of  o*/c  as  a  function  of  aa  for  Gages  615  and  625 
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5.2.3  Relationship  Between  ag  and  a/d 

Based  on  the  apparent  correlation  between  the  parameters  of  the  non-Gausaan 
probability  density  function,  aa  will  be  used  to  represent  the  general  trends  of  the  other 
two  parameters  \idien  plotted  as  a  function  of  a/d.  This  relationship  is  provided  in  Figure 
5.1 1  to  demonstrate  general  trends  between  the  wave  aa,  and  the  local  environment 
represented  in  terms  of  sea  severity/water  depth  in  non-dim^asional  form.  The  lower 
values  of  a/d  correspond  to  the  deeper  gage  locations  and  the  trend  is  for  these  values  to 
approach  the  values  defined  by  wave  profiles  having  a  normal  distribution.  This  plot  also 

r 

reinforces  the  analysis  done  in-  Chapter  4  vv^ch  determined  the  boundary  values  of  the 
Gaussian  distribution.  At  a/d  =  0.06,  the  aa  values  appear  to  begin  converging  towards 
the  Gaussian  limit  located  at  the  origin  of  Figure  5. 1 1. 

The  higher  values  of  a/d  represent  the  shallower  water  depths.  It  appears  that  at 
these  values  aa  tends  to  approach  a  finite  asymptotic  limit.  This  is  consistent  with  the 
findings  fi'om  the  spatial  analysis  figures.  In  between  these  two  limits  of  Figure  5.11  there 
is  a  depression  in  the  severe  sea  condition  case.  This  phenomenon  must  be  analyzed  more 
thoroughly  to  ascertain  a  reasonable  explanation  for  its  existence  and  will  be  the  subject  of 
a  later  study. 


0.18 
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Figure  5.11  Plot  of  aa  as  a  Fuui 


CHAPTER  6 

SUMMARY  AND  CONCLUSIONS 


It  is  shown,  based  on  the  results  of  this  study,  that  wind  generated  waves  exhibit 
significant  changes  in  their  profile  forms  as  they  propagate  into  finite  water  depths.  The 
actual  nature  of  this  transformation  has  been  a  subject  that  has  not  been  covered  in  any 
significant  detail  The  unique  wave  data  records  collected  during  the  ARSLOE  Project  of 
1980  enabled  a  detailed  study  on  this  topic  to  he  produced. 

One  objective  of  this  study  was  to  verify  the  applicability  of  the  non-Gaussian 
probabihty  density  fimction  developed  by  Ochi  and  Ahn  (1994b)  defined  as 
— exp  -  — — i — — (1  -  ya/n,  -  exp(-yax))^  -  yax 

,  L28^r:x>0 
where,  y  =  1 

^  3.00for:x<0 

Through  statistical  analysis  of  the  ARSLOE  data,  it  has  been  confirmed  that  this 
probabihty  density  fimction  can  be  used  with  confidence  to  represent  wave  profiles  in 
finite  water  depths.  It  correlates  very  well  with  histogram  data  over  an  extensive  range  of 
sea  conditions.  The  distribution  of  the  wave  profile  in  finite  water  depths  ^ows  a 
dependence  on  water  depth  and  sea  severity.  In  particular,  as  water  depth  decreases,  the 
distribution  becomes  increasingly  non-Gaussian.  Also,  as  sea  severity  increases  the  same 
trend  occurs. 

With  the  appUcabihty  of  the  probabihty  density  fimction  verified,  an  analysis  of  the 
boimdary  condition  where  the  wave  profile  can  no  longer  be  considered  Gaussian  was 
conducted.  A  general  relationship  has  been  proven  to  exist  between  the  largest  significant 
wave  height  and  the  local  water  depth  where  the  wave  profile  can  stih  be  considered 
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H  (7 

Gaussian.  This  relationship  is  defined  as  =  0.24  •  Both  forms  of  this 

relationship  are  identical  since  the  significant  wave  height,  Hs  is  defined  as  four  times  the 
standard  deviation,  a,  in  deep  water  conditions.  This  relationshq)  has  been  verified  using 
three  different  approaches  with  each  yielding  similar  results.  In  shallower  water 
conditions,  this  boundary,  defining  the  transition  fi-om  Gaussian  to  non-Gaussian  wave 
profiles,  appears  to  shows  a  dependence  on  beach  slope.  Hiis  is  a  qualitative  observation 
fi-om  the  present  study  and  additional  research  must  be  done  to  quantify  this  finding.  This 
dependence  on  beach  slope  may  be  similar  to  the  influence  of  beach  slope  on  breaking 
waves.  The  linear  result  defined  by  this  study  is  akin  to  the  earliest  breaking  criterion 
based  on  Solitary  Wave  Theory. 

At  present,  the  only  way  this  probability  density  fimction  can  be  defined  is  through 
knowledge  of  statistical  properties  of  the  wave  field  of  interest.  In  particular,  the  mean, 
variance  and  third  moment  of  the  wave  profile  must  be  known.  This  requires  that  wave 
profile  data  be  obtained  in  order  to  evaluate  the  parameters  of  the  probability  density 
fimction. 

Another  objective  of  this  study  was  to  determine  relationships  linking  the 
parameters  of  the  probability  density  fimction  to  local  field  conditions.  If  the  wave  profile 
statistics  can  be  defined  in  terms  of  local  environmental  conditions  such  as  sea  state  and 
water  depth,  the  probability  density  function  can  be  defined  without  wave  data.  This  study 
has  revealed  qualitative  trends  between  the  parameters  of  the  probability  density  function 
and  local  field  conditions  based  on  the  extensive  data  complied  fi-om  ARSLOE.  From 
these  qualitative  trends,  it  is  evident  that  the  specific  relationships  between  the  parameters 
and  environmental  conditions  are  complex  and  may  require  consideration  of  additional 
parameters  besides  sea  severity,  a,  and  water  depth.  The  pursuit  of  a  functional 
relationship  between  the  parameters  and  the  local  environment  continues  to  be  a  topic  of 
considerable  interest  and  active  research. 


APPENDIX 

COMPARISON  OF  ORIGINAL  AND  SMOOTHED  “a”  VALUES 


It  can  be  seen  that  the  nonlinearity  parameter,  “a”  is  characterized  by  a  significant 

amount  of  local  flucuations  superimposed  on  the  overall  trend  of  the  parmeter.  Based  on 

the  definition  of  “a”  as  a  nonlinearity  parameter  used  to  describe  the  probability  density 

> 

fimction  of  a  wave  profile,  these  '^nall  flucuations  can  be  smoothed  without  any  noticable 
loss  of  accmacy.  It  is  highly  unlikely  that  the  probability  density  fimction  of  a  wave  profile 
changes  so  sporadically  over  such  short  periods  of  time  as  indicated  by  the  minor 
variations  of  “a”.  There  was  not  any  one  smoothing  technique  that  was  used  in  the 
smoothing  of  “a”  at  aU  gage  positions.  Instead,  several  smoothing  techniques  were 
evaluated  and  the  technique  which  produced  the  lowest  mean  absolute  deviation  (MAD) 
was  the  technique  used  for  that  particular  gage  position. 
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